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ABSTRACT 

We observed a small sample of 5 radio-quiet QSOs with integral field spectroscopy to search for possible extended emission in the 
Lyo- line. We subtracted the QSO point sources using a simple PSF self-calibration technique that takes advantage of the simultaneous 
availability of spatial and spectral information. In 4 of the 5 objects we find no significant traces of extended Lya emission beyond 
the contribution of the QSO nuclei itself, while in UM 247 there is evidence for a weak and spatially quite compact excess in the Lyo- 
line at several kpc outside the nucleus. For all objects in our sample we estimated detection limits for extended, smoothly distributed 
Lyo- emission by adding fake nebulosities into the datacubes and trying to recover them after PSF subtraction. Our observations are 
consistent with other studies showing that giant Lya nebulae such as those found recently around some quasars are very rare. Lya 
fuzz around typical radio-quiet QSOs is fainter, less extended and is therefore much harder to detect. The faintness of these structures 
is consistent with the idea that radio-quiet QSOs typically reside in dark matter haloes of modest masses. 
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. Introduction 


The circum-galactic medium (CGM), defined as gas outside of 
the main stellar body of galaxies but still within the virial radii 
of their dark-matter haloes, is of crucial importance in galaxy 
evolution. It may act as a reservoir for fuelling star formation in 
the galaxy, and it is also subject to feedback processes that ex- 
pell material from the galaxy. If violent enough, this feedback 
may in turn heat up the CGM and thus prevent it from contribut¬ 
ing further to the formation of stars. Observing the CGM at high 
redshifts will hence provide relevant insights about galaxy for¬ 
mation. 

One observational approach to study the ‘cold’ gas phase 
of the CGM (T $ 10‘' K) at high redshift uses absorption 
signatures against background sources. This has provided im¬ 
portant statistical constraints on several properties of the CGM 


(e.g. Adelberger et al.|2005t [Hennawi et al.|2006} Steidel et al. 

2010| l. However, the spatial distribution of the CGM in individ¬ 
ual galaxy cannot be captured with this method. 

An alternative approach is to map the CGM in the Lyman 
a emission line. Several mechanisms have been identified that 
should lead to Lya emission from the CGM: Cooling of infalling 


gravitationally heated gas (e.g., Haiman et al. 200(11, cooling 
following superwind-driven shocks (e.g., Taniguchi & Shio^ 
2000)1 and - possibly most important for our investigation - Lya 


fluorescence induced by exposure to UV radiation. While the 
metagalactic UV background alone is predicted to produce only 
a very faint glow that is probably out of reach for the current gen¬ 
eration of optical instruments (e.g., Kollmeier et al.|2010|l, Lya 


* Based on observations collected at the Centro Astronomico His- 
pano Aleman (CAHA) at Calar Alto, operated jointly by the Max- 
Planck Institut fur Astronomie and the Instituto de Astrofisica de An- 
dalucfa (CSIC) 


fluorescence caused by the much stronger UV radiation from a 


QSO should boost the emission into the detectable regime (Rees 
|1988[ [Haiman & Rees|200T| |Kollmeier et al!]|2010| l. Searching 
for Lya signatures of the CGM around luminous QSOs is the 
topic of the present study. 


IHaiman & Rees) P001| l estimated that a z ~ 3 QSO har- 
bouring a 5 x IO^^Mq halo should be surrounded by Lya fuzz 
extending radially outwards to 2"-3", at a surface brightness of 
~ 5 X 10“'^ erg s^'cm^^arcsec'^. At least in this model, Lya 
fuzz is predicted as a generic property of high-z QSOs, where 
the surface brightness of this fuzz depends only on the mass of 
the halo. More recent theoretical work suggests, however, that 
a substantial fraction of gas is accreted within filamentary cold 
streams (Dekel et al.|2009 Faucher-Giguere & Keres 2011 Ros 


|dahl & Blaizot|2012[ and references therein). If these streams are 
optically thick to ionizing radiation, they will develop a highly 
ionized skin in the presence of a quasar. This skin then acts like 
a mirror converting up to two thirds of the incident ionizing ra¬ 
diation into Lya photons. In this scenario the expected surface 
brightness then depends on the ionizing photon flux produced by 
the quasar and the projected spatial configuration of the streams 
( [Kollmeier et al.|2010[|Hennawi & Prochaska|2013| ). 


Revealing extended Lya structures around QSOs requires 
a proper subtraction of the PSF-broadened nuclear component, 
which will outshine - even under good seeing conditions - the 
expected CGM signal close to the QSO. Observationally, this 
makes the detection of circum-QSO Lya fuzz much harder than 
searching for Lya ‘blobs’, now routinely found in large-area 


narrow-band surveys (Steidel et aL||2000l Matsuda et aL| 

2004 

Saito et al.||2006l Ouchi et al. 20091 IMatsuda et al.||201 

lllErbl 

et al. 1)20111 Prescott et al.||2012 

2013)1. Many of these ‘blobs’ 


have no obvious central source of ionizing photons (e.g., |Nils-| 
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|son et aL]|2006| l, which might make them physically distinct 
from circum-QSO Lya fuzz. On the other hand, some Lyo blobs 
may also be powered by AGN that are highly obscured along the 


datacubes to unveil possible extended Lya nebulae. We detail 
how we estimate surface brightness upper limits in Sect. 3.3 and 


line of sight (e.g. 

Basu-Zych & Scharf|2004HGeach et al.|2007l 

Hayes et al. 201 

1 Martin et al. 2014b|l; such blobs would be 


physically linked to the fuzz we want to observe (see also Baek 
|& Ferrara|2013| l. 

As of yet, sizeable samples of extended Lya emission around 
QSOs exist only for radio-loud objects ( [Heckman et al.|1991a|b| l, 
for which such emission seems to be a generic property at z ~ 2- 
3, with detections reported up to z ~ 6 ( [Roche et al.j 201^ . 
However, it appears likely that a large fraction of extended Lya 
emission around radio loud quasars is driven by interaction with 
the radio jets. Evidence for this stems from an observed spa¬ 
tial correlation between radio and Lya morphology (e.g., [Heck^ 
man et al.|1991b| Humphrey et al.[[2007[l and a frequent occur¬ 


rence of similar morphologies in extended C iv and He ii emis¬ 
sion ( [Humphrey et al.|2006[ [Sanchez & Humphrey|2009| l. 

In order to search for the CGM in QSOs unaffected by ra¬ 
dio jets, one obviously has to resort to radio-quiet objects. For 
this class, however, mostly single object discoveries have been 


published so far ([Bergeron et al.|1999( Bunker et al.|2003| Wei- 

dinger et a 

1 2004 2005 [ Francis & McDonnell 

2006 Willott 

et al. 

2011 

Goto et al.[ 2012 Rauch et al.[ 20 

3| Cantalupo 

et al. 

2014 

Martin et al.[[2014a|l, and only one 

single-object 


study reported a non-detection at a faint surface-brightness de¬ 
tection limit ( Francis & Bland-Hawthorn]2004[ l. Very few pro¬ 
grammes aimed at observing (however small) sam ples ([Chris- [ 
[tensen et al.[[2006| [Courbin et al.pOOS) [North et al.[[2012| l, and 

these reported high detection frequencies of circum-QSO Lya- 


fuzz: [Christensen et al.[ ( 

2006|l found Lya fuzz around 4 of 6 

radio-quiet QSOs, while 

North et al.l (|2012]l, by extending the 

sample of Courbin et al. 
6 targets. So the literatui 
rather generic feature anc 

(200811, found Lya fuzz around 4 of 
■e suggests that Lya fuzz is indeed a 
not a peculiarity around high-redshift 


present the results of those calculations. We summarise our con¬ 
clusions in Sect.[4| 

Conversions of observed to physical quantities always as¬ 
sume standard ACDM cosmology with Qa = 0.7, Qm = 0.3 
and Ho = 70kms 'Mpc"'. 


2. Observations & Data Reduction 

2.1. Sample selection & observations 

Our sample was selected to be at redshifts around z ~ 2.3 so 
that, in the case of a successful detection of extended Lya fuzz, 
spectroscopic follow-up observations of rest-frame optical emis¬ 
sion lines would be possible. Starting from the [Veron-Cetty &[ 
[Veron ( [2010| l catalog with this constraint, we browsed through 
all bright (my > 19 as in Veron-Cetty & Veron ( 2010[ l) radio¬ 
quiet QSOs at approximatly the above redshift. In the end we 
were able to observe 5 of these targets with the Potsdam Multi- 
Aperture Spectrophotometer (PMAS; [Roth et al.|2005 i at Calar 
Alto Observatory during 3 consecutive cloudless dark- to grey¬ 
time nights in October 2011. Table lists the redshifts and the 
SDSS u and g band magnitudes of our observed QSOs. 

We used PMAS in its 16x16 lens array (LArr) configura¬ 
tion, giving an 8"x8" field of view (FoV), (i.e. at 0.5"x0.5" spa¬ 
tial sampling). Motivated by the increased sensitivity and larger 
spectral coverage of the upgraded PMAS detector ( [Roth et aL 
[2010[l, we mounted the VI200 grating to obtain the highest pos- 


radio-quiet QSOs, even if the number of systematically observed 
QSOs appears still too low for certainty. 

A major increase in sample size was recently achieved due 
to the work of [Hennawi & Prochaska[p013| l. In a long-slit cam¬ 
paign on 8m-class telescopes they observed 29 close quasar- 
quasar pair sightlines, where Lyman-limit absorption in the 
background quasar spectrum indicated the presence of optically 
thick hydrogen clouds in the CGM of the foreground z ~ 2 
radio-quiet quasars. The a priori known presence of such gas 
clouds implies that the sample would be positively biased for 
detecting ‘mirrored’ Lya fuzz converted from the quasars ion¬ 
izing photons (see above). They find a single quasar surrounded 
by a large scale nebula and 10 objects surrounded by small scale 
Lya fuzz. Surprisingly, none of their spectra showed Lya fuzz 
with the properties expected for optically thick ‘Lya mirrors’. 

For a survey on circum-QSO Lya fuzz, integral-field spec¬ 
troscopy is an optimal observational method since it has con¬ 
tinuous spatial and spectral coverage. This technique allows for 
an optimal subtraction of the PSF-broadened nuclear emission. 
Moreover, since spatial and spectral information are obtained 
simultaneously in case of a detection, more inferences on the 
physical state of the CGM gas can be made. Given the previous 
success of a Calar-Alto 3.5m PMAS IFU campaign where 4 ex¬ 
tended Lya nebulae were detected around 6 radio-quiet QSOs 
( [Christensen et al.|2006| l (hereafter CJW06), we initiated a new 
targeted PMAS IFU campaign to extend this sample. 

The outline of this paper is as follows: In Sect.j^we describe 
our PMAS observations and how we reduced the raw data. In 
Sect. 1^ we explain how we remove the QSO emission from the 


sible spectral resolution within the targeted wavelength range. 
The CCD was read out unbinned in dispersion direction. Expo¬ 
sure times, seeing, airmasses and observing conditions of our 
observations are listed in Table To ensure the best possible 
spectral tracing and wavelength calibration we flanked each on- 
target exposure by continuum and HgNe lamp exposures. At the 
beginning and end of each night we observed an [Oke ( 1990| l 
spectrophotometric standard star (BDH-24d4655 and G191B2B, 
respectively). We obtained skyflats during twilight and several 
bias frames during the night in idle time (e.g. while performing 
acquisitions). 

2.2. Data reduction with p3d 

To reduce the observations we employed the p3d-pipelin^ 
( [Sandin et al.|2010 [2012| l. We briefly outline in this section how 
we applied the different tasks of the pipeline to our data. 

For every observing night a masterbias was created with 
the routine p3d_cmbias. The p3d_ldinask task then produced 
dispersion solutions for every arc-lamp frame, fitting 5* order 
polynomials for the mapping from pixel- to wavelength space; 
measured residuals for all arc lines were < 10 'px. To deter¬ 
mine spectrum traces and cross-dispersion profiles we applied 
the p3d_ctrace method to every continuum lamp frame. We 
created wavelength calibrated flatfields using p3d_cflatf. The 
p3d_cobjex task then performed extraction, flat-fielding, and 
subsequent wavelength calibration of every target and standard 
star exposure. As detailed by [Sandin et al.[ ( [2010[ l, the best 
signal-to-noise for PMAS LArr spectra is achieved by using 
the modified optimal extraction (MOX) algorithm ( [Horne[1986[ l. 
Due to the large separation of LArr fiber traces on the CCD, 
cross-talk correction between fibers is not needed. I Sandin et al.l 
P012|l advises correcting for small spectra trace shifts between 


* http://p3d.sourceforge.org 
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Table 1. Observed radio-quiet QSOs: Properties, exposure Times & observing conditions. 


QSO 

Zvc 

Zsdss 

msdss 

.gSDSS 

^exp [s] 

Seeing ["] 

Airmass 

Conditions 

Q0027-H0103 

2.310 

2.313 

19.27 

18.50 

7x1800s 

0.9-1.3 

1.2-1.6 

photometric / grey 

Q0308-H0129 

2.335 

- 

19.21 

18.45 

7x1800 

1.0-1.4 

1.2-1.5 

mainly clear / dark 

Q0256-0003 

2.381 

2.385 

20.13 

19.41 

8X1800-H 1x900 

0.9-1.1 

1.3-1.6 

photometric / dark 

Q2243-H0141 

2.314 

- 

18.48 

17.98 

8x1800 

0.9-1.2 

1.3-1.5 

photometric / grey 

UM247 

2.333 

2.333 

18.85 

18.10 

5x1800s 

0.9-1.1 

1.2-1.4 

photometric / dark 


Notes. Redshifts from |Veron-Cetty & Veron]j2010^ (zvc); and |Paris et al.|p012^ (zsdss)- Photometry from |Ahn et al.|pbl2^ . 


continuum-lamp and on-target exposures for achieving highest 
fidelity with MOX extraction, therfore we switched on median 
recentering in p3d_cobjex. Compared to extractions where we 
experimentally turned off this feature, we saw a clear improve¬ 
ment in signal-to-noise, even though the determined offset was 
typically less than half a pixel. Cosmic-ray hit removal was im¬ 
plemented within p3d_cobjex, utilizing the PyCosmic ( |Huse-| 
|mann et al.|2012[ ) algorithm. We also chose the option to subtract 
a scattered light model before extraction, as advised for PMAS 
LArr data ( |Sandin et al.|2012] l. Next, the p3d_sensfunc routine 
was utilized to create sensitivity functions. Here we first created 
a ID standard star spectra with the p3d_sens£unc GUI using 
the extracted standard star observations. Then these sensitivity 
functions were applied to every extracted target exposure using 
the p3d_fluxcal task. For correction of atmospheric extinction 
we used values from the empirical Calar Alto extinction curve 
model ( [Sanchez et al.|20()7| l. 

The final data products resulting from the p3d pipeline are 
flux calibrated datacubes for all target exposures. p3d also pro¬ 
duces the corresponding error cubes, containing the standard- 
deviation of each volume pixel (‘voxel’). After trimming regions 
affected by vignetting on the detector (Roth et al. 20101, the 
cubes cover a wavelength range from 3600-4600 A sampled 
on identical wavelength grids (Ad = 0.75 Apx^*). Their spec¬ 
tral resolution, determined from fitting ID Gaussian profiles to 
several lines in the spectra of the extracted arc-lamp frames, is 
vfwhm ~ 160 kms^' (R « 1850). We note there is a spatial 
and spectral dependence of R in PMAS (see also [Sanchez et al.] 


2012 Sect. 6.2), and the value reported here is the median near 


the centrum of our wavelength range. 

2.3. Sky subtraction 

To remove night sky emission from the datacubes, we created for 
each target a median spectrum from the ring of spaxels bracket¬ 
ing the FoV (ignoring low-transmission fibres) and subtracted 
this from every spaxel. Even though the median is more robust 
against contamination of signal from the target than the mean, 
it is still possible that we subtracted a fraction of the nebular 
emission if this emission extended out into and beyond the FoV. 
However, we note that most of the known extended Lyor region^ 
around RQQs have projected maximum extents of < 65 kpc, cor¬ 
responding to < 8" (the PMAS FoV) at the redshifts of our ob¬ 
jects. Nevertheless, we first visually inspected all datacubes if 
any obvious extended emission features were present, and then 
checked the individual outer-ring median spectra if they contain 
spikes not attributable to sky lines. Both tests were negative, so 
we are certain at this point not having accidentally removed very 


bright extended nebular emission. Nevertheless, we will return 
to this point in our analysis in Sect. [3.5[ showing that for special 
cases of large-scale extended emission our observing strategy 
may have been less than optimal. 

2.4. Stacking the individuai exposures for each QSO 

To increase the signal-to-noise of faint emission we stacked the 
individual sky-subtracted & flux-calibrated datacubes. During 
the observations the QSO was centered in the LArr FoV using 
the acquisition and guiding (A&G) system of PMAS. No dither¬ 
ing scheme was intended, to avoid cross-correlation of neigh¬ 
bouring spaxels on the sky. Unfortunately the filter wheel of 
the A&G camera was set to V-Band, causing the QSOs centroid 
point source to move around within our bluer wavelength range 
due to differential atmospheric refraction, even though the guide 
star observed with the A&G system remained steady. We cor¬ 
rected for these unintended dither offsets by spatial integer pixel 
shifts before stacking. Integer pixel shifts were preferred over 
fractional shifts to avoid interpolation effects and thus keeping 
the original observational information as unaffected as possible. 

In practice, we first determined the QSOs centroid position 
in the individual exposures using images created from summing 
up datacube layers around /iLyff(zQso)- If a centroid was shifted 
more than a half pixel along a spatial axis, this cube was then 
shifted by a pixel along this axis in the opposite direction. No 
exposure needed to be shifted by more than one pixel in x or 
y direction. Because of the varying sky background (and differ¬ 
ent exposure times in the case of Q0256-0003) we employed a 
variance-weighted mean for stacking, using the error cubes gen¬ 
erated by the pipeline. 

Having reduced the raw data as described, we were now left 
with five datacubes containing voxels where (x, y) are the 
spatial indices, and z is the index for the spectral layers, in units 
of 10“^^ erg s^^cm^^A. 


3. Analysis & results 

3.1. QSO spectra 

We show spectra of our five QSOs in Fig. [T] These spectra were 
extracted in a 3" diameter aperture centered on the highest S/N 
spectral pixel (spaxel). Where available, we overlay SDSS DR9 
( |Ahn et ar[2012| l spectra. Our spectra differ to SDSS by < 


10 

for Q0027-H0103 & Q0308-H0129; for UM247 there is a sys¬ 
tematically offset by < 20 % in the blue. Quasar variability, spec- 
trophotometric uncertainties in the SDSS or imperfections in our 
flux calibration could be reasons for these differences. We note 


^ Using the compilation in Fig. 3 of 

Cantalupo et al. 

2014 

and 10 that for Q0308-I-0129 the observing conditions were not photo- 

additional detections reported by Hennawi & Prochaska 

2013 

L metric, thus our fluxes of this object are likely somewhat too low. 
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'^Obs. [A '^Obs. [A 


Fig. 1. Extracted spectra from the PMAS datacubes in a 3"diameter aperture. For those objects where available (Q0027+0103, Q0256-0003 & 
UM 247) SDSS DR9 spectra have been over-plotted in grey to illustrate the quality of our flux calibration. The vertical dashed line shows the 
wavelength of Lya at the QSOs redshift and the vertical dotted lines indicate the wavelengths of the artificial telluric Hg i (4047 A & 4358 A) 
emission-lines. 


The night sky emission line spectrum in Calar Alto shows 
significant man-made contributions, arising from tropospheric 
scattering of high-pressure street lamps in nearby populated ar¬ 
eas. In our wavelength range the Hg 1 4047 A & 4358 A emission 
lines are prominent. Wavelengths of these lines and the expected 
Lya wavelength 4°^® = (1 +ZQSo)'^Lyii indicated in Fig. As 
can be seen, the 44047 line is located unfavorably close near the 
expected peak of the Lycr emission for 4 of our QSOs. We had 
not anticipated the strength of this anthropic line in the planning 
of our observations. This unresolved line has an averag e surface 
brightness of = 3.5 x 10“'® ergs“*cm“^arcsec“^ (Sanchez 
|et al.|2007] l and thus, by amplifying the background noise in a 
narrow-band window around 4°®^-, contributes negatively to our 
efforts in uncovering extended Lya emission. The sensitivity in 
this narrow-band is further decreased by residuals from sky sub¬ 
traction, resulting from the varying spectral resolution across the 
FoV (Sect.f 


3.2. Subtraction of OSO emission from the datacubes 

Due to the PSF broadening, spectra at radial distances < 2" from 
the QSO center are visibly contaminated by the nuclear spec¬ 
trum. Thus, in order to reveal possible extended Lya emission 
around the QSO, we had to remove this contamination. |CJW06| 
experimented with three algorithms to achieve this deblending 
in PMAS LArr datacubes; 

- On- minus scaled off-band method: An on- and off-hand im¬ 
age is created from the datacube. The latter is then scaled 
and subtracted from the former, thereby revealing possible 
contribution of extended emission in the on-band. Here the 
on-band is obtained as a summation in wavelength direction 
over all datacube layers expected to contain the extended- 
emission signal, while the off-band is a sum over nearby lay¬ 
ers not being affected by extended emission. Asserting that 
the extended emission contributes only a few percent to the 
on-bands peak intensity (i.e. QSO-i-extended emission), the 
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ratio between the peak intensities in on- and off-hand is a 
sufficient scale factor for the subtraction. 

Analytical PSF Extraction: An analytical PSF model is fitted 
and subtracted from every datacube layer, thereby creating 
a datacube which should only contain non-nuclear emission. 
Then an iterative scheme is employed, where in a second it¬ 
eration positional and shape parameters of the PSF function 
are fixed, using values from the first iteration with the con¬ 
straint of allowing them to vary only smoothly with wave¬ 
length (see e.g. Wisotzki et al.|2003) [Kamann et al. 2013 


for more elaborate versions of this technique involving mul¬ 


tiple point-sources within the FOV). According to CJW06 


the analytical PSF extraction produces unsatisfying results 
in PMAS in terms of QSO residuals, so we did not consider 
applying this method to our data. 

Empirical PSF Subtraction: By summation over datacube 
layers expected to be unaffected by extended emission, an 
empirical PSF image is created. It is then scaled and sub¬ 
tracted from every layer, thus creating a residual datacube 
containing extended emission. To find the scaling the PSF 
image with respect to the datacube layers, a spectrum is 
aperture-extracted from the datacube, with every spaxel in 
the aperture having weights assigned via the PSF image. An 
iterative scheme is employed to remove contamination in the 
scaling spectrum from extended emission possibly covering 
parts of the aperture: the datacube resulting from the de¬ 
blending is subtracted from the undeblended datacube from 
which in turn the next scaling spectrum is extracted. A vari¬ 
ant of this method was developed by |Husemann et ak ( 2013) 1, 
where in subsequent iterations the correction is achieved us¬ 
ing an annulus-extracted spectrum from the deblended dat¬ 
acube. 


We emphasize the methodological similarity between the on- 
minus scaled off-band method and the empirical PSF subtrac¬ 
tion deblending: If the PSF image in the latter is created from 
the same layers as the off-hand in the former, every layer in the 
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deblended datacube can be thought of as an ‘on-band’ having 
a scaled ‘off-band’ subtracted. If, moreover, no iterative correc¬ 
tions are applied, the image resulting from the on- minus scaled 
off-band method is identical to the image resulting from summa¬ 
tion over the on-band defining layers in the empirical PSF sub¬ 
tracted datacube. Note that the iterative correction only increases 
the fidelity of an extended emission signal when it is measurably 
present after the first deblending iteration; otherwise only noise 
is shuffled around in subsequent iterations. 

Because of this similarity we applied solely the empirical 
PSF subtraction method to our datacubes. As PSF image we take 
~ 40 layers (~ 30 A) beginning ~10 A redwards of dLyaCzQso)- 
We found this choice to produce the best results in terms of sub¬ 
traction residuals, because of the following reasons: Firstly, for 
the red side of /iLya(ZQSo), a redward band gives a better signal- 
to-noise ratio for the PSF image due to Zabs ~ Zem Lya absorbers 
appearing blueward of /iLyff(zQso) (see Fig. [^. Secondly, as we 
did not correct for differential atmospheric refraction, the posi¬ 
tion and shape of the PSF changes with A, thus selecting layers 
further away would produce a non-optimal representation of the 
PSF at the expected position of Lya fuzz, resulting in stronger 
residuals from the subtraction (dLyaCzQso))- This is also the rea¬ 
son why we chose the spectral window to be relatively narrow. 

We experimented with different extraction apertures for the 
scaling spectrum. Visual inspection revealed no extended emis¬ 
sion after the first iteration in all objects. As expected (see 
above), performing more iterations does not improve upon this. 
We thus decided to use the smallest possible aperture for scaling, 
i.e. the single spaxel with the highest S/N. Despite this, we could 
not find any evidence for extended Lya fuzz around any of our 
targets, except a small-scale feature to the north of UM 247. 

To visualize our findings, we show in Fig. narrow-band 
images centered on /lLya(2Qso) created by summing over 20 lay¬ 
ers (15A) of the residual datacubes. Instead of using physical 
units, the scale on the colour bar uses multiples of the standard- 
deviation per pixel. This colour bar scaling simplifies the judg¬ 
ment whether features seen in this image are actually significant, 
and it also makes the comparison between the panels straight¬ 
forward. Except for UM247, none of these images shows any 
significant feature. 

In Fig. we show spectra from the residual datacubes ex¬ 
tracted within a small circular aperture consisting of 20 spaxels 
(aperture with outer radius rc - 1.25", cf. Sect. |3.3| centered 
around the scaling spectrum. Again, only UM 247 displays a dis¬ 
tinct spectral line that seems inconsistent with noise. We discuss 
this feature further in Sect. 14.21 


3.3. Estimation of detection iimits 

From our visual inspection of the residual datacubes (see also 
Figs.|^&|^ we conclude that in 4 out of 5 objects we have no 
significant detections. We now want to constrain an upper limit 
in surface brightness for extended Lya nebulae around those 
QSOs. To do so, we define confidence limits for rejecting the 
null hypothesis, which states that there is no extended Lya emis¬ 
sion present in the data. The upper limit is then given as the 
minimum surface brightness for which we would not be able to 
confidently reject the null hypothesis anymore. 

We adopted circular apertures centered on the scaling spec¬ 
trum. In absence of morphological information for the non¬ 
detections, this simplification appears reasonable, although it is 
known that QSO CGM Lya nebulae can be asymmetric (e.g. 
[Weidinger et al.|2004[ [Rauch et al.|2013| l. For such nebulae the 


average surface brightness within a circular aperture will gen¬ 
erally be lower than a surface brightness obtained within an 
isophote. 

We explored successively larger apertures around the central 
scaling spectrum by adding annuli with a width of one pixel. 
This defined 8 “circular” apertures Q (k = 1... 8) with outer 
radii of = (2k H- l)/4". The numbers of spaxels Nc in those 
apertures are then Vq = (8, 20, 36, 68, 96, 136, 176). 

The assessment was performed on the residual datacube 
after PSF subtraction. The residual signal Sf, within the 
relevant spectral layers and inside each circular aperture Ck can 
be written as 

(x,y)£Ct 


where - Yiz&i'Rp'xyz value of a pixel in the narrow- 

band images shown in Fig. For the sake of brevity, we do not 
explicitly state units and conversions in our equations. 

Assuming for the moment that we had reliable variance esti¬ 
mates cr^,, for every pixel of the pseudo narrow-band image 
(or equivalently cr\y^ for every voxel of the residual datacube), 
the noise in Ct could be written as 



( 2 ) 


Combining now Eqs. ([T]) and Q, we could express the signal Su 
within Ck in amounts of the noise CTk being present in Ck, i e. 


Sk ^n- a-k , 


(3) 


where, under assumption of pure Gaussian noise, n would di¬ 
rectly translate into a probability of the null hypothesis - no neb¬ 
ular emission detected - being false (e.g. |Wall|1979^ . 

However, neither of the above made assumptions - hav¬ 
ing reliable variance estimates for every voxel and pure Gaus¬ 
sian noise - are met by our data: The QSO subtraction intro¬ 
duces non-Gaussian residuals, as does the background subtrac¬ 
tion (Sect. [273] l, which are not captured by the formal variances. 
There may also be further unknown systematics. We therefore 
replaced the formal variances by two empirical proxies for <Tx,y. 

- The standard deviation (relative to an expectation value 
of zero) per pixel of the narrow-band image Ix,y outside of the 
circular aperture Ck. 

- The standard deviation (relative to an expectation 
value of zero) per spaxel within the circular aperture Ck, cal¬ 
culated from spectral layers not contributing to the narrow- 
band images. The calculation is limited to wavelengths not 
further than 75A away from /I(zqso), since here the empirical 
PSE subtraction starts to produce strong residuals because of 
differential atmospheric refraction. 

Note that both noise estimators depend on the the aperture Q. 

With these proxies we then set ctx.j ~ or ~ 17^'’“ 
for all X, y in Eq. Q and obtain 

Sk = «spaf(i?r/V^)’ (4) 

Sk = «spec ■ (<“/#q) . (5) 

We show below that Uspec ~ Wspat (or equivalently 17^'’^”^) 

holds. 
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Q2243+0141 UM247 



Fig. 2. Residual spectra for our 5 QSOs from the QSO-subtracted datacubes, extracted within the = 1.25" aperture (cf. Sect, around the 
scaling spectrum. The vertical dotted lines indicate the boundaries of the narrow-band image (15 A = 20 layers), shown in Fig.l^ The vertical 
dashed-dotted lines show the position of the Hg i-sky line (Sect. ED- 
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Fig. 3. Narrow-band images (15A wide, centered on /iLyaizQso)) for our 5 QSOs. These images were created from the QSO-subtracted datacubes 
using the empirical PSF subtraction method. The position of the spaxel used for scaling the PSF image is indicated by a cross; this position 
corresponds to the centroid of the QSO nucleus. Colours indicate multiples of the standard-deviation per pixel CTpi^, calculated for each image. The 
values of are (2.3, 1.4, 2.7, 2.6, 3.5)x 10^'^erg s“' cm“^, from top left to bottom right. North is up, and East is to the left. Axes ticks are AS and 
Aa in arcseconds with respect to the QSO centroid. 


The question is, -which rtspec or Hspat is required, in each QSO, 
for a detection? We addressed this problem by adding simu¬ 
lated extended emission into our datacubes before subtracting 
the QSO. Specifically -we used circular nebulae of an extent that 
would fill a particular aperture Q. These simulated nebulae had 


a flat surface brightness profile and a Gaussian line profile with 
300 km s ' FWHM (approximately twice the spectral resolution) 
centered around /lLyQ.(ZQso)- We emulated seeing effects by con¬ 
volving those nebulae with 2D Gaussians of the average seeing 
FWHM of the particular observation. Using the above defined 
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Fig. 4. Examples of recovered spectra after adding artificial nebulae into the datacube prior to the empirical PSF subtraction. The surface brightness 
of the artificial nebulae was scaled with integer values of Uspatiai according to Eq. After empirical PSF subtraction the spectrum was extracted 
from the residual datacube using the rc - 1.25" aperture (i.e. the same as in Fig.® The spectra are shown in units of 10^'^erg s^'cm^^A^'. Note 
that artificial nebulae with Uspa, > 5 can be unambiguously discriminated from the background for all objects. 



Fig. 5. Comparison of formal 5 <t surface brightness detection limits as a 
function of aperture radius using the different noise estimators ITspat and 
(Tspec, exemplarily shown for two objects (Q0027-I-0103 - blue symbols, 
UM 247 - red symbols panel). The agreement between ITspat and ITspec is 
similarly good for the other objects. 


noise proxies we scaled the surface brightness of the nebulae by 
integer multiples Wspat according to Eq. ffl. We downsampled 
our simulated nebulae to the grid of our aatacubes and added 
them before the final step of empirical PSF subtraction. By vi¬ 
sual inspection of the residual cubes we found that nebulae with 
Hspat = 5 can be unambiguously discriminated from the back¬ 
ground. Exemplarily we show in Fig.|^the results of this numer¬ 
ical experiment for nebulae covering the = 1.25" aperture. 
Results for other aperture sizes were similar, i.e. a input 

according to Eq. Q yielded an unambiguous visual detection 
after PSF subtraction. 

We now demonstrate that both noise proxies yield similar re¬ 
sults for our surface brightness limits, defined as the minimum 
surface brightness that a circular nebulae with a particular radius 
could have before it would fall below our detection criterion. 


For this purpose we show in Fig. the surface brightness lim¬ 
its exemplarily for two objects, as a function of aperture radius. 
These limits were calculated using either Eq. Q or Eq. (|^ with 
«spat = 5 or Hspec = 5. Figure shows the trend expected for a 
noise estimate independent of aperture, obtained by scaling the 
curve for the smallest aperture (rc = 0.75") by aperture area 
oc Note that indeed and have similar values, 

and that there is only a mild dependence on aperture size. For 
the cases not shown, the agreement is similar. Since also no a 
priori distinction can be made as to which of the two proxies is 
better, we quote the average -H (T'’’’“‘)/2 as our detection 

limits in surface brightness. 


3.4. Surface brightness limits on extended emission 


The resulting surface brightness upper limits for extended Lya 
emission surrounding the observed QSOs, calculated by apply¬ 
ing the method presented in the previous section, are shown in 
Fig. 0 For reference, the obtained upper limits for Lya fuzz 
with 2.5" (20.5 kpc) radial extend - the typical extent of circum- 
QSO Lya fuzz predicted by |Haiman & Reei] ( |2001) l - are (5.4, 
3.4, 6.5, 6.0, 6.3) x 10“^^ erg s“^ cm“^ arcsec ^ for the quasars 
Q0027-H0103, Q0256-0003, Q0308 h-0129, Q2243-h 0141, and 
UM 247, respectively. 


In Fig. 1^ we also show the integrated signal with a circu¬ 
lar aperture, i.e. Eq. Q transformed to surface brightness units. 
Here the error bar on the integrated signal indicates the standard 
deviation within the aperture. In our case, this is a measure for 
how asymmetric possible signal is distributed within the aper¬ 
ture; e.g. for UM247, a bright spot appears only to the north of 
the object, thus the error bar on the integrated signal within the 
circular aperture is large. Note that for all objects except UM 247 
the integrated signal is always below our detection limits, thus 


confirming the visual impression gained from Figs. 


discuss the fuzz around UM 247 further in Sect. 14.2 


^andig We 


Background offsets, positive or negative, can be seen in the 
integrated signal for Q0308-I-0129 or Q2243-I-0141 - however, 
note that we incorporated these systematics in the calculation of 
our detection significances by forcing the expectation value to 
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Fig. 6. Surface brightness upper limits (triangles) for circularly symmetric extended Lya emission around our 5 QSOs as a function of radial 
angular (bottom abscissa) or physical extent (top abscissa). The integrated observed signal within each apertures is shown by the black points, 
where the error bars indicates the standard deviation within the aperture, as a measure of the irregularity of the flux distribution within the aperture. 
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Fig. 7. Expected average surface brightness profiles for exponen¬ 
tial nebulae SBLypfr) = Eo x exp(-r/ro), measured within circular 
apertures and compared to the our surface brightness upper limits 
{dashed lines, indicating the range of our limits shown in Fig. [^. 
Coloured solid lines show the expectations including the effect of 
our sky subtraction procedure, while the dashed lines ignore it {blue: 
Eo ~ 10^'^ergs^* cm^^arcsec^^ and = 2"; green: Eq = 5 x 
10^'^ergs^' cm^^arcsec^^ and ro = 4"). 


zero for our noise proxies (see above). There also appears a hint 
of possible extended emission in Q0256-0003, which is how¬ 
ever below our detection threshold and not confidently separable 
from noise in the residual datacube. 


3.5. Effect of sky subtraction on large-scale Lya fuzz. 


Recently, two nebulae were discovered that have sizes signifi¬ 
cantly larger than our field of view ( Cantalupo et al.|2014[ Mar^ 
tin et al. |2014a see Table [^. As explained in Sect. ll' we 


subtracted a median spectrum generated from the outer spax- 
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els framing our field of view to remove emission by the night 
sky. Would giant nebulae as the above mentioned still leave a 
detectable signal in our observations? 


From the intensity map presented in Cantalupo et al. ( 2014| l 
we see that at a distance of r ^ 5" to the QSO UM 287, their neb¬ 


ula has a mean surface brightness of a; 5 x 10 
arcsec 


erg_s 


cm 


At a similar distance to their central QSO, Martin et al. 

9 X 1 


erg s 


( 2014a I measured a surface brightness of : 
cm“^ arcsec^^. If the surface brightness profile increases towards 
the center, the central parts might still leave a significant de¬ 
tectable signal in our datacubes. If, however, the central surface 
brightness profile is rather fiat, the significance of the recov¬ 
erable signal will be substantially reduced. Unfortunately, cur¬ 
rently no information on the central parts of the surface bright¬ 
ness profile in UM 287 is available (Cantalupo, priv. comm.). 


We investigated the recoverability of the nebula in UM 287 
assuming two different exponential surface brightness pro¬ 
files, SBLya'(r) = 5^0 X exp(-r/ro) with ro = 2" and 

ro = 4", which are fixed to UM287’s surface brightness 
of lO^'^erg s“'cm^^arcsec^^ at r = 5". The central sur¬ 
face brightness of the ro = 2" profile is then Eo ~ 

10“'®ergs“* cm^^arcsec^^ and for the ro = 4" profile Eq = 5 x 
10^'^erg s“^ cm^^arcsec^^. Our sky subtraction procedure would 
subtract a constant surface brightness of SB(4") a; 1.5 x 10“'^ 
from both profiles. We note that such a faint signal would not 
have been seen in the visual inspection of the subtracted sky 
spectra. 


Applying our detection criterion in Fig.]^ the ro = 2" profile 
would still permit a significant detection of the nebula. However, 
for the flatter ro = 4" profile the recovered signal would fall be¬ 
low the detection threshold. And obviously we would always un¬ 
derestimate the true extent of the nebula. We emphasize that our 
observations were originally not designed with such large scale 
nebulae in mind and their recent discovery came as a surprise. 
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Table 2. Compilation of observational results from the literature on extended Lya emission around radio-quiet QSOs (single-object investigations). 


Reference 





Quasar 

Redshift 

Extent 

Flux 

Method 

Bergeron et al.|( 

19991 


Q2233-606 

2.238 

9.2"xl2.1" 

3.2x10-^^ 

narrow-band imaging 

Bunker et al. 

(2C 




PC0953-H4749 

4.46 

=s5" 

3.6x10-'^ 

long-slit spectrum 

Weidinger et al. 

(2005 


Q 1205-30 

3.041 

=sl0" 

7x10-1® 

narrow-band & long-slit 

Francis & McDonnell 

(2006 

1 PSS 21554-1358 

4.28 

<0.7" 

1.7x10-1^ 

integral field spectroscopy 

Goto et al.|(2 

UT^ 

.1 



CFHQS J2329-0301 

6.417 

-1.3 -1.5" 

8x10-1^ 

long-slit spectrum 

Rauch et al. 

201 

3i 



QJ0332-2751 

3.045 

2.2" 

5.4 X 10-1^ 

long-slit spectrum 

Cantalupo et al. 

2014: 

UM287 

2.28 

^55" 

5.5x10-1® 

narrow-band imaging 

Martin et al.]f20 

14a 

1 


HS 14494-19 

2.843 

>30" 

2.0x10-1® 

integral-field spectroscopy 


Notes. Flux is given in erg s“*cm“^. Measurements on the fuzz around CFHQS J2329-0301 have also been published by Willott et al. l j2011 1 , who 
obtained similar results. No flux has been reported by |Bunker et al.| ( |2003[ l for the flux of the extended Lya emission surrounding PC 0953-1-4749, 
but this object is also part of the |CJWO^ sample. The flux for the nebula around HS 1449-1-19 is from Table 3 in |Martin et al.| ( [2^4a^ . See text for 
possible caveats considering the comparability of the presented values on flux and extent. 


Table 3. Compilation of observational results from the literature on extended Lya emission around radio-quiet QSOs (samples). 


Reference Radio-Quiet Quasars Redshift SB limit Instrument - Method 


Targeted Fuzz detected 


CJW06 l"l 





6 

4 

2.7 < z < 4.5 

- 10-11 

PMAS - integral-field spectroscopy 

North et al. 

2012 

, (b) 



6 

4 

z~4.5 

- 5 X 10-1® 

FORS2 - slit spectroscopy 

Hennawi& 

Prochaska 

(2013 

1 

29 

11 

z ~ 2 

- 5 X 10-1* 

LRIS & GMOS - slit spectroscopy 


This study 5 1 z ~ 2.3 ~ 5.5 x 10 PMAS - integral-field spectroscopy 


Notes, (a) Surface brightness limit for |CJW06| estimated by us. (b) |North et al.|p012| l is an extension of |Courbin et al.|f2008| l. Surface brightness 
limit in erg s“*cm“"arcsec“^. 


4. Discussion 

4.1. Comparison with observations from the literature 

How do our upper limits in surface brightness for extended Lya 
emission around radio-quiet QSOs compare with previous in¬ 
vestigations of this phenomenon? Most of the studies reporting 
a successful attempt in detecting circum-QSO Lya fuzz (some¬ 
times serendipitous discoveries) have focused on single objects. 
We compiled a list of such investigations in Table Beyond 
these single-object results, very few studies aimed at construct¬ 
ing actual samples (which were always small). We list the rele¬ 
vant publications in Table 

Substantial methodological differences between the studies 
listed in Tables [^ & [^ have to be kept in mind when compar¬ 
ing those results with our non-detections. The sizes in Table 
often refer to the maximum extent at which the authors were 
able to detect Lya emission. This quantity obviously depends on 
the depth of the observations and, in the case of long-slit spec¬ 
troscopy, the orientation of the slit. For the latter case, fluxes are 
also affected by significant slit-losses, since only a fraction of 
the nebula is usually captured. Finally, line emission from the 
central QSO might contaminate the fluxes of the nebular com¬ 
ponent in some cases, especially since a subtraction of the QSO 


a different setup than we) and knowing the typical instrumental 
and atmospheric parameters, we estimated that their detection 
sensitivity was similar to that in our study. 

We next considered whether the Lya fuzz detected by the 
studies listed in Tables |2] and |3] would have been recovered 
if ‘implanted’ into our quasars. For this exercise we selected 
only objects for which the reported maximum extent (after de- 
redshifting) would be covered by our field of view. The Martin| 


|et al.|(|2014a|l,|Cantalupo et al.|(|2014|l and|Bergeron et al.' 


(19991 


objects do not fulfill this criterion (see also Sect. |3.5[). We also 


point source was not performed in some cases (e.g. Bergeron 
|et al. 1 1999)[Martin et al.|2014a[ [Rauch et al.|2013| l. 

The quoted values for the surface brightness limits of the 
samples in Table are also very rough estimates, as the actual 
limits depends on the assumed size of the fuzz (see our deriva¬ 
tion in Sect. 13.31 and also the derivation in Sect. 4.3 of IHen-l 
|nawi & Prochaska|2013| l and, moreover, they differ from target 
to target due to differences in the used instruments and observ¬ 
ing strategies. Unfortunately, |CJW06| did not quantify the depth 
of their observations. Since they also used PMAS (although in 


excluded the Rauch et ar| ( |2013] l object, since their flux value is 
contaminated by quasar Lya emission. We then circularized the 
nebulae, i.e. we assumed the flux to be distributed symmetrically 
around the quasar with a radius defined as half the maximum ex¬ 
tent. After calculating de-redshifted radii and surface brightness 
levels, we determined the integrated signals within circular aper¬ 
tures covering the whole nebulae. Note that only |Weidinger et ah] 
P005| l provided a surface brightness profile, so that only for this 
object we could calculate the integrated signal at various radii. 

In Fig. we plot the results from this exercise and compare 
them to our surface brightness limits: Out of 11 nebula used in 
this calculation, 3 would be detected in all our observations, fur¬ 
ther 4 would yield a 5cr detection only in our deepest dataset, 
and 4 would not be recoverable at all. (Note however that the 
circularization is actually reducing the signal from sources that 
have significantly asymmetric flux distributions. Thus a nebula 
yielding a 5cr detection only in our deepest dataset might have 
been recovered at a higher significance if the area over which the 
signal was integrated was better matched to the light distrubution 
of the nebula.) 

Note that half of the points in Fig are at radii less than 
10 kpc and 8 out of 10 are below 25 kpc. Small scale Lya fuzz 
appears to be quiet frequent around radio-quiet QSOs. This is 
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Fig. 8. Comparison of our surface brightness upper limits as a func¬ 
tion of radius (thick dashed lines) to the reported literature detections 
of circum-QSO Lyn emission, de-redshifted to z = 2.3, and assuming 
that the reported maximum extent defines the radius of the detection 
aperture. The symbols feature results from the |CJW06| sample (blue 
triangles with error bars), from the North et ^j2012] isample (green 
squa res), the nebula from|Goto et al. ^2012^ (red diamond) and the fuzz 
from [Francis & McDonnell| ( |2006| l (cyan circle). We also show the in¬ 
tegrated signal at various radii from the de-redshifted jWeidinger et'H^ 
([2005^ surface brightness profile (solid green line). 


also supported by the high recovery rate of this phenomenon in 
[Hennawi & Prochaska| ( |2013l l, although their detections are gen¬ 
erally fainter. We argue in the next section that the signal we find 
around UM 247 also falls into this category. 

Finally, we point out that most of the reported detections of 
bright Lycr fuzz are typically at z > 3. Our recovery rate of this 
phenomenon at z ^ 2.3 (20%) is lower compared to CJW06 
(66%). Moreover, Hennawi & Prochaska ( 2013| l report a recov¬ 
ery rate of 38% at z ~ 2. Observations are insufficient at this 
stage to infer a redshift evolution of circum-QSO Lya fuzz prop¬ 
erties. It is nevertheless intriguing that Zirm et al. ( |2009 1 find 
that Lya nebulae enshrouding radio-loud QSOs decline in size 
and luminosity with decreasing redshift. If real, such a decline 
could indicate a depletion of the cool gas content of the COM 
(at least on average). Extending such a trend to radio-quiet QSOs 
would however require much larger and more homogenous sam¬ 
ples over a range of redshifts. 


4.2. Small-scale Lya fuzz near UM247 

UM 247 is the only object in our sample with a formally signifi¬ 
cant residual at Lya after PSF subtraction. To assess whether this 
residual is real and not an artefact caused by a fault (e.g. an unde¬ 
tected cosmic) in a single exposure, we repeated the stacking of 
the individual exposures (Sect. |2.4| i 5 times, each time with one 
exposure excluded from the stack. We then performed the em¬ 
pirical PSF subtraction in exactly the same manner as described 
above (Sect. [3^ for each of those stacks. The feature remained, 
hence we assert that it is genuine. 

This basically unresolved excess emission is located at a dis¬ 
tance of d <4-8 kpc (0.5-1") north of UM247. The line has a 
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flux of ~ 5.6 X 10 erg s 'em corresponding to a Lya lumi¬ 
nosity of Liya ~ 2.4 X lO''^ erg s“'. 

If the Lya emission of this feature was powered purely by 
star formation, ignoring ionization of the nearby quasar as well 
as radiative transfer effects, this would correspond to a star for¬ 
mation rate of ~ 10-20 Mq yr“' (with the exact value being de¬ 
pendent on the metallicity of the underlying stellar population, 
see Schaerer|2003| l. This Lya luminosity is the bright end of the 
Lya emitter luminosity function for z ~ 2.3 ( Blanc et al.|2011) l. 
Given the small field of view of our observations, the detection 
of such a bright Lya emitter in close vicinity to a quasar would 
appear rather coincidental. 

In an opposing scenario we consider the possibility that the 
Lya radiation is produced purely by fluorescence of a com¬ 
pact nearby cloud, devoid of internal star formation and opti¬ 
cally thick to ionizing radiation from the quasar. This cloud will 
behave as a special mirror, converting ~66% of all impinging 
hydrogen-ionizing photons into Lya. Thus, given the quasar ion¬ 
izing luminosity and the size of the cloud one can predict its Lya 
luminosity (e.g. [ Francis & Bland-Hawthorn ||2004 Adelberger 
jet al.|200^|Kollmeier et al.|2010[]Hennawi & Proc iaska|2013| lT 


Since the cloud is unresolved, we adopt as an upper limit of its 
size the extent of one PMAS spaxel. Thus the physical surface 
extent of the Lya emitting cloud is ^16 kpc^. From a mea¬ 

sured QSO flux density of - 2.5 x 10“'® erg s^'em^^A^' for 
UM 247 at 4500A, we estimate a quasar luminosity at the Lyman 
edge of = 8.0 x 10^® erg s“'Hz“', assuming a power-law in¬ 
dex of ai = 0.44 (fy oc v““) for the quasar continuum redward 
of the Lyman edge ([Vanden Berk et ar||2001|l, and 02 - 1.57 


blueward of 912A ( Telfer et al. 2002|l . We also corrected for 


galactic extinction (A^ - 0.08 at 4500A for UM 247), although 
that is a small effect. Assuming isotropic radiation, this results 
in an ionizing photon number flux of <1) = 4 x lO'" s ' cm“^ 
(1 X 10 '® s“' cm“^) at a distance of li = 4kpc (8 kpc). The ex¬ 
pected Lya luminosity of a spherical cloud at this location fol¬ 
lows as 


^hya — X fgai X /7thick ^ ItViya X X O , (6) 

where //thick is the fraction of ionizing continuum photons con¬ 
verted to Lya photons, i.e. //thick = 0.66. With /gm we denote 
the geometric reduction factor - a free parameter that accounts 
for the inhomogeneous illumination of the cloud and subsequent 
redistribution of Lya photons over a wide solid angle. Radiative 
transfer simulations suggest that for a transversely illuminated 
cloud /gm = 0.5 (|Kollmeier et al.|2010|). For our upper limit on 
the cloud radius. Equation (|^ then provides an upper limit on its 
Lya luminosity of Liya ^ 4 x 10"'^ erg s“' (1 x 10“'"' erg s“') at 
d = 4 kpc (8 kpc). This is almost an order of magnitude higher 
then the observed value. In reality, however, the cloud size might 
be much smaller than our instrumentally imposed upper limit, 
and it might also be further away than the projected transverse 
distance. Hence this order of magnitude estimate of Lya radia¬ 
tion emanating from the surface of an optically thick cloud in the 
vicinity of the quasar is still consistent with what we observe. 


4.3. Comparison with models 

Haiman & Rees| ( [2001[ ) presented a strongly idealized model to 
predict the luminosity of Lya fuzz around radio-quiet quasars. 
Specifically, they assumed a spherical symmetric 2-phase gas 
distribution in pressure equilibrium within a collapsed dark- 
matter halo and predicted extended Lya emission as a generic 
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property of high-z QSOs. The phases are a hot tenuous virial- 
ized plasma (i.e. Thot ~ Tyir of the halo) and colder neutral gas 
that has cooled down to Tcoid ~ 10"^ K during the age of the 
system. If in such an environment the quasar emits radiation 
isotropically and ionizes the whole nebula, its Lya luminosity 
depends only on the total gas mass and thus on the total halo 
mass. In this framework, the absence of significant Lycr fuzz in 
our QSOs might suggest that the hosting haloes are not overly 


massive. However, as already pointed out by Haiman & Rees 
( |2001| l and also |Alam & Miralda-Escude| ( |2002[ ), small devia- 
tions from this idealized scenario might alter the surface bright¬ 
ness of the Lya fuzz substantially. Currently, the observations do 
not provide strong constraints on the basic assumptions in these 
models. 

State-of-the-art numerical simulations in a cosmological 
framework predict that the spatial distribution of the CGM 
gas shows filamentary structure, with cold gas accreting along 
streams towards the center of the halo (e.g. Dekel et al. 2009 [ 


|Faucher-Giguere & Keres|201 1 |Rosdahl & Blaizot 2012| l. Those 

streams have high column densities and their surfaces will there¬ 
fore reflect up to 2/3 of incoming ionizing photons as Lya pho¬ 
tons. Although Lya cooling radiation from those streams alone 
might already produce a detectable signal in extremely massive 
haloes, the presence of a central quasar should enhance the con¬ 
trast of the filamentary structures by boosting their Lya emis- 
sivity by up to 2 orders of magnitude ( Kollmeier et al.||2010| l. 
While the giant Lya nebulae around some quasars may well be 
explained by such fluorescently glowing accretion streams, these 
are by no means typical. More high-resolution simulations will 
be required in order to predict robustly the luminosities and sizes 
of Lya fuzz for the lower-mass haloes typical for radio-quiet 
QSOs. 


5. Conclusions 

It seems that on average, radio-quiet QSOs are rather unspec¬ 
tacular sources of spatially extended Lya emission. The nonde¬ 
tection of such Lya fuzz in 4 of our objects and the marginal 
detection in one case are all fully consistent with the results of 
other recent investigations, although even the combined samples 
are still small. 

The spectacularly bright and extended Lya nebulae discov¬ 
ered ar ound a few quasars (|Weidinger et al.|2005[ [Martin et al.| 
|2014a[ [Cantalupo et al.|2014| l must be probably considered very 
rare cases. The rarity of this phenomenon may be explained if 
giant Lya nebulae are seen only around QSOs that reside in ex¬ 
traordinarily massive haloes. 

The observing techniques used to search for Lya fuzz around 
QSOs are quite divers, encompassing narrow-band imaging, 
long-slit spectroscopy, and now also integral held spectrocopy 
(IPS). In principle, IPS should surpass the other methods by a 
large margin; in fact, any IPS datacube allows the user to ex¬ 
plore both the narrow-band imaging as well as the spectral do¬ 
main. The limiting factor for most existing optical IPS studies 
- including the present investigation - is sensitivity and light¬ 
collecting area of the available telescopes. This is however about 
to change, as a new generation of efficient IPS systems is being 
deployed at 8-lOm class telescopes. Of particular interest for 
the topic of this study is the MUSE instrument, recently com¬ 
missioned at the ESO Very Large Telescope ( Bacon et al.|2014| l. 
Its unprecedented sensitivity will make it an optimal discovery 
machine for Lya fuzz around quasars. Moreover, the large field 
of view of MUSE will ensure tfiat tfiere is no longer a danger 
of sky subtraction removing physical signal from very extended 


nebulae with flat radial profiles. It is to be expected that within a 
relatively short time, the statistics of observed Lya fuzz around 
quasars will improve dramatically, turning the emphasis from 
discovery to the detailed dissection of physical properties. 
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